Non-adiabatic effects in long-pulse mixed-field orientation of a linear polar molecule 
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We present a theoretical study of the impact of an electrostatic field combined with non-resonant 
linearly polarized laser pulses on the rotational dynamics of linear molecules. Within the rigid rotor 
approximation, we solve the time-dependent Schrodinger equation for several field configurations. 
Using the OCS molecule as prototype, the field-dressed dynamics is analyzed in detail for experimen- 
tally accessible static field strengths and laser pulses. Results for directional cosines are presented 
and compared to the predictions of the adiabatic theory. We demonstrate that for prototypical 
field configuration used in current mixed-field orientation experiments, the molecular field dynamics 
is, in general, non-adiabatic, being mandatory a time-dependent description of these systems. We 
investigate several field regimes identifying the sources of non-adiabatic effects, and provide the field 
parameters under which the adiabatic dynamics would be achieved. 

PACS numbers: 37.10.Vz, 33.80.-b, 33.57.-|-c,42.50.Hz 



I. INTRODUCTION 

The ereation of direetional states of moleeules repre- 
sents an important tool to control and tailor the rota- 
tional degree of freedom. When a molecule is oriented 
the molecular fixed axes are confined along laboratory 
fixed axes and its dipole moment is pointing in a particu- 
lar direction. Experimentally, the availability of oriented 
molecules provides a wealth of interesting applications 
in a variety of molecular sciences, such as in chemical 
reaction dynamics , photoelectron angular distribu- 
tions , or high-order harmonic generation [§, [13] ■ 

Due to this broad interest, special efforts have been 
undertaken to create samples of oriented molecules and 
techniques based in the application of inhomogeneous [11, 
[TTI | , and homogeneous [l3 - [l6l e lectric fields as well as ho- 
mogeneous magnetic fields |l7l| have been used. A major 
breakthrough came with the proposal by Friedrich and 
Herschbach [3 [3 of enhancing the orientation of polar 
molecules by exposing them to combined weak electro- 
static and strong non-resonant radiative fields. This the- 
oretical prediction was done within an adiabatic picture 
assuming that the switching on time of the laser pulse 
is larger than the molecular rotational period. For lin- 
ear molecules, a linearly polarized laser field produces a 
double- well potential along the polarization direction. In 
the pendular limit, this double-well potential contains 
nearly degenerate pairs of states with opposite parity 
forming tunneling doublets. If the molecules possess a 
permanent electric dipole moment, a strong pseudo-first- 
order Stark effect is induced by coupling the tunneling 
doublets with an additional electrostatic field. Due to 
this coupling, the two levels in a pendular doublet are 
eSicicntly oriented but with their effective electric dipole 



moments pointing in opposite directions. As a conse- 
quence of this oriented and antioriented states pairing, 
the orientation is small in a molecular ensemble with the 
population thermally distributed. Therefore, the first ex- 
perimental measures of the orientation of a molecular 
beam were indeed reduced to small values [l^, [2l|. A 
significant improvement was gained by using a quantum- 
state selected molecular beam, which allowed the cre- 
ation of unprecedented degree of orientation for com- 
plex asymmetric tops [22| - f24j . A first theoretical study 
of the mixed-field orientation experiment of asymmetric 
top molecules, pointed out that a fully adiabatic descrip- 
tion of the process does not reproduce the experimental 
observations [25j . 

Recently, we have experimentally and theoretically in- 
vestigated the mixed-field orientation of the carbonyl sul- 
fide molecule (26j . Our analysis has proven that a time- 
dependent description of the mixed-field orientation pro- 
cess is required to explain the experimental results. We 
have shown how the non-adiabatic coupling of the levels 
forming the quasi-degenerate doublets as the laser in- 
tensity is increased, gives rise to the reduction of the 
orientation and, therefore, to the disagreement with the 
predictions of the adiabatic theory [3, [3 ■ Herein, we 
provide a detailed theoretical analysis on the dynam- 
ics of a linear molecule exposed to an electrostatic field 
combined with a non-resonant laser pulse. In the frame- 
work of the rigid rotor approximation, we solve the time- 
dependent Schrodinger equation using experimental field 
configurations, i.e., a Gaussian laser pulse and a weak 
electrostatic field that is turned on at constant speed. 
As prototype example, we consider the carbonyl sulfide 
molecule (OCS). For several rotational states, we investi- 
gate the mixed-field orientation dynamics under different 
field-configurations by varying either the laser peak in- 
tensity, the duration of the Gaussian pulse, the dc field 
strength or the angle between both fields. Hence, we 
demonstrate that for some field configurations, the field- 
dressed dynamics is non-adiabatic and provide a detailed 
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account of the sources of non-adiabaticity and the field 
regimes at which they appear. For parahel fields, the 
dynamics is characterized by the population transfer be- 
tween adiabatic states when the pendular doublets are 
formed. Whereas for non-parallel fields, we encounter 
additional non-adiabatic effects when the states from the 
same J-manifold, having now the same symmetry, are 
driven apart as the laser intensity is increased on the 
weak field regime. For different field configurations, we 
identify and discuss the experimental conditions needed 
to achieve an adiabatic molecular dynamics. 

The paper is organized as follows: In ISec. Ill we de- 
scribe the Hamiltonian of the system and its symmetries 
for various field configurations. The results for the en- 
ergy, alignment and orientation predicted by the adia- 
batic theory arc analyzed in lSec. lUl In lSec. IVl we focus 
on the molecular dynamics when the fields are parallel. 
In particular, we explore how the time-dependent ori- 
entation varies as the field parameters are modified, and 
indicate the experimental conditions under which an adi- 
abatic orientation would be achieved. A similar study is 
performed for tilted fields in ISec. VI where we show that 
the conditions for an adiabatic mixed-field orientation are 
more difficult to fulfill. In ISec. VH we assume that once 
the pulse is turned on its peak intensity is kept constant, 
and investigate the dynamics in this regime. In ISec. VIII 
first the laser pulse is switched on and then an electric 
field is applied. In this field configuration, we analyze 
the orientation of the ground state and provide the field 
parameters for an adiabatic orientation. The conclusions 
are given in ISec. Villi 

II. THE HAMILTONIAN OF A LINEAR RIGID 
ROTOR IN EXTERNAL FIELDS 

We consider a polar linear molecule exposed to an ho- 
mogeneous static electric field and a non-resonant lin- 
early polarized laser pulse. The field configuration is 
illustrated in |Fig. 1| the polarization of the laser field 
lies along the Z-axis of the laboratory fixed frame (LFF) 
{X,Y,Z), and the dc field is contained in the XZ-plane 
forming an angle /3 with the .Z-axis. The z-axis of the 
molecule fixed frame {x^y, z) is defined by the permanent 
dipole moment fi of the molecule. These two frames are 
related by the Euler angles fl = {0,(j)), cf |Fig. 1| The 
description of this system is done within the rigid rotor 
approximation, assuming that the vibrational and elec- 
tronic dynamics are not affected by the fields. Thus, the 
rigid rotor Hamiltonian reads 

H{t)=Hr + H,{t)+Hi^{t), (1) 

where Hr is the field-free Hamiltonian 

Hr = B3\ (2) 

with J being the total angular momentum operator and 
B the rotational constant. The terms Hs{t) and HL{t) 



stands for the interactions with the static and laser fields, 
respectively. 




FIG. 1. Laboratory fixed coordinate, Euler angles, schematic 
field configuration and the OCS molecule. 

The dipole coupling with the static field reads 

i/,(i)=-/^-E,(i) = -^E,(t)cos0, (3) 

with Es(t) = Es(t)(sin/3X-t-cos^Z), and Es(t) being the 
electrostatic field strength. The angle between the dipole 
moment /x and this field is Bs , cf |Fig. 1[ and cos Og = 
cos /3 cos + sin j3 sin cos 0. 

We consider a non-resonant laser field linearly po- 
larized along the Z-axis of the LFF, EL(i) = 
Eo5(t) cos(27ri'i)Z, with v being its frequency, Eq the 
peak field strength, and g{t) the pulse envelope. Assum- 
ing that is much shorter than the pulse duration and 
the rotational period, we average over the rapid oscilla- 
tions of the non-resonant field. This causes the coupling 
of this field with the permanent dipole moment to van- 
ish [23, [13. Thus, the non-resonant laser field molecule 
interaction can be written as 

ifL(<) - -^Aacos^fl, (4) 
2ceo 

where Aa is the polarizability anisotropy, is the in- 
tensity of the laser, c is the speed of light and eo is the 
dielectric constant. Note that we have neglected the term 
— a_Ll(t)/2ceo, which represents only a shift in the energy. 

In this work, the field configurations are chosen based 
on the mixed- field orientation experiments [13, [13, [1^ ■ 
Initially, the molecule is in field-free space, then the elec- 
trostatic field is switched on, and its strength is increased 
linearly with time. At t = —Tq the maximum strength 
is achieved and kept constant afterwards. This time Tq 
is chosen long enough to ensure the adiabaticity of this 
turning-on process. For the laser pulse, we use a linearly 
polarized Gaussian pulse with a full width half maximum 
(FWHM) T on the nanosecond range. The intensity is 
given by l{t) = Iq exp (— i^/2cr^), with Iq being the peak 
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Table I. Action of the symmetry operations on the Euler angles. 



Transformations 



Operation (/> 9 

1 4^~T^ e~Te 

Cx{i^) (j>-^2-K-(j} 9^ TV -9 

Cz{S) <j>^<f) + 5 9^9 



Table II. For the /3 = 90° and 0° < < 90° field configurations, 
values of the parameters k and t for the wave functions with the 
correct symmetry for each irreducible representation. 
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intensity, and a is related with the FWHM t = 2\/2 hi 2(7. 
Numericahy, the non-resonant laser field is turned on in 
such a way that the interaction due to this field is much 
weaker than coupling with the dc field. 

The eigenstates of the field-free Hamiltonian |(2)| are 
the spherical harmonics Yjm{^), with J and M being 
the rotational and magnetic quantum numbers, respec- 
tively. Note that M is the projection of the total angular 
momentum J on the LFF Z-axis. The field-free Hamil- 
tonian (2) belongs to the SO (3) group because the op- 
erator remains unaltered under any rotation. In the 
presence of the external fields, the symmetries of the ro- 
tational Hamiltonian (1) arc significantly reduced. The 
operations describing the symmetries of the field-dressed 
Hamiltonian are collected in ITable II For a static field, 
the symmetry group is reduced to arbitrary rotations 
around the field axis and the identity {£ ,Ces(^)}- If 
only a linearly polarized laser field is applied, the Hamil- 
tonian is invariant under arbitrary rotations around Z- 
axis and two-fold rotations around any axis perpendicu- 
lar to the Z-axis, and the symmetry group is compound 
by {£ ,C z{5) tC _\_z{t^)} ■ For parallel fields, the elements 
of the symmetry group are the identity £, arbitrary ro- 
tations around the fields Cz{S) and the reflection in any 
plane containing the flclds. For a given |Af|-value, the 
parity under the reflection on one of these planes deflnes 
two irreducible representations. Since the selection of 
this plane is not unique, then the states with M 7^ 
are doubly degenerated. Note that M remains as a good 
quantum number for these three field configurations with 
/3 = 0°. For non-perpendicular and non-collinear fields, 
i.e., /3 ^ 0° and f3 ^ 90°, the Hamiltonian is invari- 
ant under the identity E and the reflection on the XZ- 
plane containing the fields axz- This group {£,(7xz} 
has only two irreducible representations characterized by 
the parity with respect to this reflection axzi and the 
functions belonging to the even and odd representations 
arc ^y, ^(ri) = { YjMm + {-irYj.Min)) /V2, with k 
given in ITable III If the fields are perpendicular, /? = 90°, 
the two-fold rotation around the static field Cx (tt) is also 
a symmetry operation. Thus, the symmetry group is 
formed by {£, cxZiCjc (tt)} has four irreducible represen- 
tations labeled by the parity with respect to the trans- 
formations axz and Cx{Tr)- The properly symmetrized 
functions for these four irreducible representations are 
^^Y^^in) = {YjMin) + {-iyYj_Mm)/V2, the pos- 
sible values of e are collected in ITable Til 



The time-dependent Schrodingcr equation associated 
to the Hamiltonian (1) is solved by means of a second- 
order split-operator technique (29j . combined with the 
discrete- variable and finite-basis representation methods 
for the angular coordinates [30|-|33|. For reasons of ad- 
dressability, we will label the time-dependent states as 
\J,Mjf and |J,M,Z) for /3 = 0° and 0° < /3 < 90°, 
respectively, with I = e and o indicating even or odd 
parity with respect to the XZ-p\anc. The labels J and 
M refer to field-free quantum number to which they are 
adiabatically connected. Note that the labeling of the 
states depends on the way the fields are turned on [33 |. 
The time-dependent wave function depends on the time 
t, the peak intensity Iq, the FWHM r, the electrostatic 
field strength Eg, and the angle (3. For the sake of simplic- 
ity, we have not made explicit these dependences, but the 
field configuration is clearly indicated through the text. 

To get a better physical insight on the field-dressed 
dynamics, the time-dependent results will be compared 
to those from the adiabatic theory. For this system, we 
take the adiabatic limit by using a constant electrostatic 
field Eg and constant laser intensity I in the Hamiltonian 



(1) The time- independent Schrodinger equation associ- 
ated to this Hamiltonian is solved by expanding the wave 
function in a basis that respects the symmetries. The adi- 
abatic states are labeled as \J,MJ)° and |J,M,Op for 
/3 = 0° and 0° < /3 < 90°, respectively, and we have not 
made explicit their dependence on the field parameters. 

The field-dressed eigenfunctions of this time- 
independent Hamiltonian form a basis, which is 
used to analyze the time-dependent wave function I7) 
by means of the following expansion 



N 



i7) = E^'Wi^' 



(5) 



with Ci{t) — (7|7i)p, 7 and 7^ including all the labels 
identifying these levels. For computational reasons, we 
have only considered the lowest-lying adiabatic levels, 
and always ensured that the contributions of highly ex- 
cited states are negligible. Let us remark that for each 
time t, the expansion of the wave function is performed in 
a different adiabatic basis obtained by solving the time- 
independent Schrodinger equation using the static held 
strength and laser intensity at time t, i.e., Es(i) and I(i). 
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III. RESULTS IN THE ADIABATIC LIMIT 

In this work, we use the OCS molecule (see fig. [T]) as 
benchmark to illustrate our results. The rotational con- 
stant of OCS is _B = 0.209 cm~^, the permanent dipole 
moment /i = 0.709 D and the polarizability anisotropy 
Aa = 4.04 A3. 



0.75 




0.25 i 



2 5 7.5 

I (10" W/cm^) 

FIG. 2. (Color online) Adiabatic results for the (a) energy, ex- 
pectation values (b) (cos^ d) and (c) (cos 6) as a function I of 
the adiabatic states |0,0, e)^ (red thick solid), |l,0,e)p (gold 
thin solid), |l,l,e)p (orange thick short-dashed), |2, l,e)p 
(dark blue long-dashed), |2,0, e)p (blue dot-short-dashed), 
|3,0,e)p (purple dotted), [2,2, e)p (green thin short-dashed), 

and |3, 2, e)p (pink dot-long-dashed). The insets show the rel- 
evant energy and intensity ranges where the formation of the 
near-degenerate doublets occurs, (d) Polar plots of the square 



of their wave functions at I 
and 13 = 0° for all data. 



lO'-'W/cm^ Es = 300V/cm 



We start by analyzing the adiabatic limit. We re- 
strict this study to the following eight states: |0, 0, e)° 



|l,0,e)°, |l,l,e)°, |2,l,e)°, |2,0,e)°, |3,0,e)°, |2,2,e)°, 

and |3,2,e)p. For f3 ~Q°, they adiabatically correspond 
to the states forming the four first doublets. Note that 
they well represent the main physical features observed 
in the overall molecular dynamics, and similar behavior 
and properties are, therefore, obtained for states in other 
irreducible representations. 

For Eg = 300V/cm and /3 = 0°, the energies and the 
expectation values (cos^ 6) and (cos 9) of these levels are 
plotted versus the laser intensity in |Fig. 2| (a), (b) and 
(c), respectively. The weak static field breaks the field- 
free degeneracy in the magnetic quantum number and. 



as the laser intensity is increased, these states become 
high field seekers. In the strong laser field regime, once 
the pendular regime is reached, pairs of quasi-degenerate 
states with the same symmetry are formed. The insets 
in this panel show how these doublets appear. The gap 
in energy in a doublet goes as AE « 2 j/iEg p(i| cos6'|i)p|, 
where fj. p{i\cos9\i)p is the effective dipole moment of 
the state |z)p in the doublet, which is of opposite sign for 
\j)p- Within a doublet, the two levels are characterized 
by the same hybridization of the angular motion (J^) and 
alignment (cos^B), see |Fig. 2| (b). In contrast, they pos- 
sess opposite orientation (cos 9) , one being oriented and 
the other antioriented, cf. |Fig. 2| (c ). This opposite ori- 



entation is also illustrated in Fig. 2| (d) by the polar plots 
of the square of their wavefunctions for I = 10^^ W/cm^. 
The larger is the field-free rotational quantum number 
of the levels, i.e., their field-free energy, the stronger is 
the laser intensity needed to achieve a significant orienta- 
tion. Indeed, the states in the third and fourth doublets 
are not aligned for I < 2 x lO^^W/cm^ and, therefore, 
not oriented. Once the pendular regime is achieved, the 
orientation of these states | (cos 9) \ approaches to 1 as I is 
enhanced. If the laser field is sufficiently strong, this adi- 
abatic orientation is independent of the dc field strength, 
and of the angle between both fields. 



IV. RESULTS FOR PARALLEL FIELDS 

In this section, we investigate the rotational dynamics 
in a parallel configuration: a dc- field of 300 V/cm and a 
Gaussian pulse with FWHM t = 10 ns and several peak 
intensities. For the ground state |0, 0,e)'^, the expecta- 
tion value (cos 9) is presented in |Fig. 3| (a) as a function 
of I(i) up until the peak intensity Iq is reached. For com- 
parison, the adiabatic results are also shown. 

Since the FWHM is 125.31 times larger than the ro- 
tational period, one would expect that the rotational 
dynamics follows the adiabatic limit. However this is 
not the case, and there are significant discrepancies be- 
tween the time-dependent and adiabatic results. In con- 
trast to what is predicted by the adiabatic theory, the 
final orientation decreases as the peak intensity of the 
laser pulse is increased. For Iq = 2xlO^^W/cm^, 
{cos 9) initially resembles the adiabatic behavior, but 
it achieves a maximum value (cos0) = 0.899 for 
l(t) = 1.86 X 10" W/cm^. For Iq = 5 x 10" W/cm^, 
10^^ W/cm^, and 2 x 10^^ W/cm^, the orientation shows 
a qualitatively similar but quantitatively different be- 
havior: in the weak laser field regime, (cos 9) mono- 
tonically increases following the adiabatic limit, but for 
l{t) > 2 X 10^^ W/cm^ it reaches a plateau behavior 
being the orientation smaller than the adiabatic value. 
For instance, (cos 6*) = 0.661 for 1(0) = 2 x lO^^^/cm^ 
whereas the adiabatic value is {cos9) = 0.975. 

A first physical insight into the non-adiabatic dynam- 
ics could be gained by analyzing the characteristic times 
of the molecule. When the states in a pendular dou- 
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FIG. 3. (Color online) (a) For the ground state, time 
evolution of the expectation value {cos 8) as a function of 
l{t) for Gaussian pulses of r = 10 ns and peak intensities 
lo = 2 X 10" W/cm^ (red thick solid), 5 x lO" W/cm^ (or- 
ange dashed), 10^^ W/cm^ (gold dotted) and 2 x 10^^ W/cm^ 
(green dot-dashed). The adiabatic results for {cos 9) (thin 
solid) are also included, (b) The squares of the projections 
of the time-dependent wave functions onto the adiabatic pen- 
dular state |0, 0, e)p. (c) Adiabatic criteria as a function of 
l{t). The field configuration is Es = 300 V/cm and ^ = 0°. 



blct arc quasi-dcgcnerated, the energy gap between them, 
AE ^ 2/iEs, defines a time scale of this system [2^1 • Note 
that wc have assumed |p(i| cos0|i)p| ~ |p(j| cos0|j)p| w 
1. which holds in the strong laser field regime. For 
Eg = 300 V/cm and (3 = 0°, the energy separation 
within the first doublet formed by |0, 0, e)p and 1 1, 0, e)p is 
AE = 6.97 x lO"** cm~^ giving a time scale of 761.21 ps, 
which is larger than the rotational period 79.8 ps. Thus, 
only long enough pulses compared to this pendular time 
would ensure an adiabatic orientation of the molecule. 

Since the static field strength is so weak, its impact 
on the rotational dynamics is very small, and before the 
pulse the levels could be considered as field-free rotor 
states. As the laser intensity increases, the states are hy- 
bridized by the combined action of the both fields, and 
the doublets of nearly-degenerate states are formed in 
the strong laser field regime, as it is shown in |Fig. 2| (a). 
When the energy splitting of this pendular doublet ap- 
proaches the coupling of the two sublcvels due to the 
pseudo-first-order Stark interaction, these states can mix 
because they have the same symmetry for /3 ^ 90° . As a 
consequence, there is a population transfer between the 
oriented and anti-oriented states, which results in a de- 
crease of the final orientation compared to the adiabatic 
limit. For this field configuration, the dynamics can be 
analyzed by means of the adiabatic states forming this 
pendular doublet, because their couplings to states in 



neighboring doublets is much smaller than the energy dif- 
ference between them. Note that these adiabatic states 
are the eigenstates of the Hamiltonian at fixed time t. 

Under a time-dependent interaction, i.e., in our case 
the interaction with the laser field Hj^{t) [(4)| the dynam- 
ics could be considered as adiabatic if and only if the 
following condition [s^ 



V = 









(■ 

p \ 


dt 


' p 



\E, - K- 



< 1 



(6) 



is fulfilled, with \i)p and |j)p being the eigenstates of the 
Hamiltonian in the adiabatic limit, and Ei and Ej their 
energies. According to this criterion, the probability for 
mixing, corresponding to a transfer of population from 
one state of the doublet to the other, is determined by the 
rate of change of the laser field interaction and the energy 
separation between the states. Thus, as the laser inten- 
sity is increased the population transfer between the two 



states in a doublet takes place because the criterion |(6) is 
not satisfied. To illustrate this phenomenon, we show the 
contribution of the adiabatic ground state |Cooe(i)P to 
the time-dependent wave function of |0, 0, e)", [Fig~3l (b), 
and the adiabatic parameter 77 when 77 > 10~^, |Fig. 3| (c). 
Note that |Cooe(i)P + 1^106(^)1^ = 1. In these four cases, 
the dynamics is initially adiabatic, i. e., |Cooe(i)P remains 
equal to 1 and rj ^ 1. As l{t) is increased, the energy 
splitting of the doublet decreases and, moreover, it be- 
comes comparable or even larger than the rate of turning- 
on the pulse; thus, 77 is close to 1, and the population 
transfer takes place. This region where rj is not negligi- 
ble corresponds to the formation of the quasi-degenerate 
doublet. Once the doublet is formed, AE = \Ei — Ej\ 
reaches a small value and slowly decreases as l{t) is en- 
hanced; but the two states are oriented in opposite direc- 
tions and their wave functions do not overlap. Therefore, 
the coupling due to the alignment laser is much smaller 
than AE, 77 ^ 1 and the population transfer does not 
take place any longer, i.e., |Cooe(OP remains constant as 
l(t) is enhanced. The larger is this population transfer, 
the smaller is the orientation compared to the adiabatic 
prediction. Since these adiabatic states contributing to 
the dynamics are quasi-degenerated and have very close 
values of the alignment and hybridization of the angular 
motion, the lack of adiabaticity is not refiected on the 
time evolution of the energy, (cos^ 9) or (J^). 

For this field configuration, the molecular dynamics 
of excited states present analogous features for (cos 6*), 
(cos^ 6) and (J^) as those encountered here for the ground 
state. 

The adiabaticity of the field-dressed dynamics is deter- 
mined by the rate of change of the laser field interaction 
compared to the largest time scale of the system. In the 
pendular regime, the energy splitting in a doublet goes 
as AE ^ 2/iEs; then, the population transfer decreases 
if Es is increased. On the other hand, by increasing the 
FWHM augments the time scale on which the pcndu- 
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lar doublets are formed, and facilitates the adaptation of 
the molecule to this field. That is, the mixed-field orien- 
tation will be more adiabatic when either longer pulses 
or stronger static electric fields are used. Let us remark 
that with the expression being the dynamics more adi- 
abatic we mean that for a certain state, the weight of 
its corresponding adiabatic state in the time-dependent 
wave function is closer to one during the dynamics. 



A. Influence of the peak intensity lo 

Here, we investigate the orientation at the maximum 
of the laser pulse, as it is done in most of the experiments 
[20I [23, ■ The rat e of change of the laser field and the 
adiabatic parameter | (6) | depend linearly on Iq. Then, for 
a Gaussian pulse with fixed FWHM, the dynamics will 
be more diabatic if Iq is increased. In this section, we 
consider the states |0,0, e)°, |I,0,e) 



|l,l,e)", |2,l,e)°. 



|2,0,e)°, |3,0,e)°, 



|2,2,e)°, and 13,2,6)". Their orienta- 



tion at t = 0, i.e., (cos 6') for 1(0) = lo, is plotted as a 
function of Iq in |Fig. 4| (a) and (b) for = 300V/cm 
and 600V/cm, respectively. The fields are parallel, and 
the FWHM of these pulses is 10 ns. 
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FIG. 4. (Color online) Expectation value {cos 6) a.t t — as 
a function of the peak intensity Iq for the states |0, 0, e)" (red 
thick solid), |1,0, e)'' ^gold thin solid) 11,1,6)" (orange thick 
short-dashed), |2, 1, e) (dark blue long-dashed) |2, 0, e)" (blue 
dot-short-dashed), 13,0,6)" (purple dotted), |2,2,e)" (green 
thin short-dashed) and |3, 2,e)" (pink dot-long-dashed), for 
^ = 0° and (a) Es = 300 V/cm and (b) 600 V/cm. The 
FWHM of the laser pulses is 10 ns. 

For Eg ~ 300 V/cm, the orientation of the low-lying 
level in a doublet increases as Iq is enhanced, reaching 
a maximum and smoothly decreasing thereafter. This 
is counterintuitive to what is expected in the adiabatic 
limit; namely, a larger orientation when the laser inten- 
sity is increased. The maximum in the orientation is 
achieved with an alignment pulse that already gives rise 
to a non-adiabatic dynamics. However, due to the cou- 
pling between the populated adiabatic states in the pen- 
dular pair p(i| cos0|j)p, the orientation is enhanced com- 
pared to what happens at the adiabatic limit. By further 
increasing Iq , the population transferred between the two 
states is enhanced, but now the coupling between them 
is very small or even zero due to their opposite orienta- 
tion. As a consequence, the final orientation decreases as 



Iq is increased. For a certain pendular doublet, the up- 
per state is antioriented, and (cos 9) shows the opposite 
behavior as a function of Iq. Regarding the third and 
fourth doublets, the states are not oriented nor aligned 
for ac pulses with Ig < 2.4 x lO^^W/cm^. Compared 
to low-lying states, their orientation is smaller and the 
maximum of (cos 0) appears at larger peak intensities. 

By increasing the static field strength to 600 V/cm the 
energy gap of the pendular pair is also increased, being 
the characteristic time scale of the system reduced. Thus, 
for the same laser pulse, the dynamics is more adiabatic, 
i.e., less population is transferred, and the final orienta- 
tion is increased, see |Fig. 4| (b). The orientation (antiori- 
entation) of the pendular states also achieves a maximum 
(minimum), but it is so shallow that it is hardly appre- 
ciated on the scale of this figure. 
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FIG. 5. (Color onhne) Projections of the time-dependent 
wave functions onto the corresponding adiabatic states as a 
function of the peak intensity Iq for the states 10,0,6)" (red 
thick soUd), |1,1,6)" (orange thick short-dashed), |2,Q, e)" 
(blue dot-short-dashed) and |2,2,6)" (green thin short- 
dashed), for (a) Es = 300 V/cm and (b) 600 V/cm. We use 
10 ns laser pulses and P — 0° . 

To illustrate the field-dressed dynamics, we plot in 
|Fig. 5| the weights of the adiabatic states associated to the 
oriented levels in these pendular doublets. For the corre- 
sponding antioriented levels, the contributions of its asso- 
ciated adiabatic state are identical to the one presented 
here, e.g., for the ground state, we present the contri- 
bution of the adiabatic ground state |Cooe(0)P, which is 



identical to the weight |Ci 



for If, 0, e). In an adia- 



batic molecular dynamics, these coefficients are equal to 
one. Note that in the considered regime only the two adi- 
abatic states of the pendular doublet contribute to the 
dynamics. For all these levels, |Ci(0)p decreases, i. e., the 
dynamics is less adiabatic, as Iq is enhanced. By increas- 
ing Eg, AE is increased; thus, under the same Gaussian 
pulse the population transfer is reduced, i.e., |Ci(0)p is 
closer to one, and the range of peak intensities that could 
be considered as adiabatic is increased. 



B. Influence of the FWHM of the laser pulse 

The duration of the Gaussian pulse plays an impor- 
tant role in the molecular dynamics. It has been shown 
that an alignment pulse of 10 ns is not enough to achieve 
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an adiabatic mixed-field orientation for molecules such 
as OCS, benzonitrile and iodobenzene 0, [H, HI, [37j . 
It has been pointed out the need of increasing the ris- 
ing time of the laser pulses to achieve the highest possi- 
ble orientation [H, [sl, By increasing the FWHM, 
the time needed to form the pendular doublets is also 
increased. For a given field configuration, at the point 
where the pulse reaches a certain intensity the adiabatic 
parameter |(6)| is reduced if r is increased. Hence, the 
molecular dynamics becomes more adiabatic and, there- 
fore, the population transfer to other pendular states is 
reduced. Here, we investigate how the directional prop- 
erties of OCS depends on the laser-pulse FWHM. For 
the same set of states as in the previous section, |Fig. 6| 
shows the orientation at t = as a function of r. The 
fields are parallel, and we consider the peak intensities 
lo = 2 X 10" W/cm2 and IQ^'^^^/cm^, and a dc field of 
Es = 300V/cm. 




2 5 10 15 20 2 5 10 15 20 

T (ns) T (ns) 



FIG. 6. (Color online) Expectation value (cos 6) at t = as 
a function of r for (a) lo = 2 x 10^^ W/cm^ and (b) lo = 
10^^ W/cm^ The fields are parallel and = 30 0V/cm. The 
states and their labels are the same as in |Fig. 4] 



The degree of orientation of the two states in a given 
pendular pair shows the same behavior as a function of r, 
but with their dipolc moment pointing in opposite direc- 
tions. In the first two doublets and Iq = 2 x 10^^ W/cm^, 
I (cos 0) I increases with t till it reaches a plateau-like be- 
havior. The second pair satisfies that |(cos0)| « 0.428 
for T > 10 ns. Since the states in the third and 
fourth doublets have not achieved the pendular regime 
for lo = 2 X 10^^ W/cm^, increasing the pulse duration 
does not have any impact on their orientation, and (cos 9) 
keeps a constant value close to zero as r is increased. For 
lo = lO^^W/cm^, the degree of orientation of all the 
states increases and approaches the adiabatic limit as r 
is enhanced. For the ground state and t = 20 ns, we ob- 
tain {cosff) = 0.913, which is very close to the adiabatic 
value (cos 9) = 0.964. 

These results show that for parallel fields, the mixed- 
field orientation dynamics of any state could be adiabatic 
if a sufficiently long pulse and sufficiently strong fields are 
used. For a 50 ns Gaussian pulse with lo = 10^^ W/cm^, 
and Eg = 300V/cm, the dynamics can be considered as 
adiabatic for aU these states, with |Ci(0)p > 0.999. 



C. Influence of the electrostatic field strength 

Since the energy splitting in a pendular doublet is pro- 
portional to the static field strength, the degree of adia- 
baticity in the molecular orientation should increase if Eg 
is enhanced, i. e., the characteristic time scale of the sys- 
tem is reduced. In |Fig. 7] we present the final orientation 
at i = of these eight states versus Eg. We have consid- 
ered two laser pulses of r = 10 ns with peak intensities 
lo = 2 X 10" W/cm2 and lO^^ W/cm^, and /3 = 0°. 
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FIG. 7. (Color online) Expectation value {cos 6) at f = as 
a function of E^ for (a) lo = 2 x lO" W/cm^ and (b) lo = 
10^^ W/cm^. We use a 10 ns laser pulse s and / 3 = 0°. The 
states and their labels are the same as in |Fig. 4\ 

For the lowest laser intensity, the orientation of the 
|0, 0, e)" and |1, 0, e)° states is constant and independent 
of the static field for E^ > 500V/cm with |(cos6')| = 
0.915. The orientation of the levels |1, 1, e)° and |2, 1, e)'' 
monotonically increases as Eg is enhanced, and we obtain 
I (cos 6*) I = 0.792 for Eg = 2kV/cm. This peak intensity 
is not large enough for the states in the third and fourth 
doublets to be in the pendular regime. Thus, these pairs 
are weakly oriented even if a strong dc field is used, e. g., 
for Eg 2kV/cm, \{cos9)\ = 0.265 and 0.153 for the 
third and fourth doublets, respectively. 

For the strong peak intensity, all the states are in the 
pendular regime, and their | (cos 9) \ increases as Eg is in- 
creased reaching a constant value for sufficiently strong 
static fields. Their orientation approaches the adiabatic 
limit and for Eg > IkV/cm, \{cos9)\ = 0.949 for the 
states of the first doublet, and 99.91% of their popula- 
tion is on the corresponding adiabatic pendular state. 
For Eg = 2kV/cm, the states in the fourth doublet sat- 
isfy I (cos 6*) I = 0.885 and 99.99% of their population is 
on the corresponding adiabatic level. 

In conclusion, by combining sufficiently strong elec- 
trostatic fields with standard Gaussian pulses, i.e., with 
experimentally accessible peak intensities of 10^^ W/cm^ 
and 10 ns FWHM, a significant orientation is obtained 
even for excited rotational levels. It is worth remarking 
that the fields have to be parallel; then, techniques such 
us the ion imaging method [l^l could not be used to mea- 
sure the degree of orientation; whereas techniques, such 
as time of flight [13, are feasible. 
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V. RESULTS FOR NON-PARALLEL FIELDS 

In this section we investigate the rotational dynamics 
when the electrostatic field forms an angle 0° < /3 < 90° 
with the linearly polarized laser pulse. The azimuthal 
symmetry is lost, and the number of irreducible repre- 
sentations is reduced to two, see ISec. Ill Thus, states 
with different field-free magnetic quantum numbers arc 
now coupled by the electrostatic field. 

The field-free wave function of the initial state is con- 
structed as an eigenstate of the operators Ces(7'') and 
axz, see ITable II i.e., | J, A'/, e)=i?y(/3)| J, M, e)°, where 
Ryiji) is the rotation operator of an angle fi around the 
LFF y-axis [i^. This ensures that this wave function 
has the correct symmetries, and that its time evolution 
corresponds, in the adiabatic limit, to an eigenstate of 
the field-dressed Hamiltonian at any time. 

Before the pulse is turned on, an important feature of 
the ground state is that its energy gap to the next state 
with the same symmetry is proportional to the rotational 
constant _£?, which is much larger than the coupling in 
the weak laser field regime. In this regime, it evolves as 
an isolated state, and its interaction to the neighboring 
levels could be considered negligible. Hence, analogously 
to the parallel field configuration, the formation of the 
doublets in the pendular regime is the only source of non- 
adiabatic effects in its field-dressed dynamics. Note that 
the lowest lying level of the odd irreducible representation 
will show the same behavior. 




(l0'2 W/cm^) I(t) (l0'2 W/cm^) 

FIG. 8. (Color online) For the ground state, evolution of the 
expectation values (a) (cosd) and (b) (cosSg) as a function 
of l{t) of a 10 US Gaussian pulse. The field configurations 
are /3 = 30° (thick) and 45° (thin) with peak intensities lo — 
2 X 10" W/cm^ (dashed), 5 x lO" W/cm^ (dot-dashed), and 
10^^ W/cm^ (dotted). The static electric field is fixed to Eg = 
300V/cm. The adiabatic results (solid) are also included. 

For the ground state, the time evolution of the expec- 
tation values {cos 9) and (cos 0s) are presented as a func- 
tion oil{t) till the peak intensity is reached in |Fig. 8| (a) 
and (b). The Gaussian pulse has 10 ns FWHM and 
peak intensities Iq = 2 x lO^^W/cm^, 5 x 10^^ W/cm^, 
and lO^^W/cm^. We consider the inclination angles 
/3 = 30° and 45°, and a dc field of = 300V/cm. 
For comparison, the adiabatic results are also included, 
with {cos 9) being independent of (3. For a certain laser 
pulse, increasing the inclination angle towards 90° im- 
plies a decrease of the energy splitting in the first dou- 
blet, AE ~ 2/iEs cos f3, and, therefore, an increase of the 



adiabatic parameter |(6)[ Note that in the energy split- 
ting we have not considered the component of the dc field 
along the LFF X-axis. Compared to the f3 — 0° config- 
uration, the dynamics could be considered as less adia- 
batic being characterized by a larger population trans- 
fer to the other adiabatic state in this pendular pair. 
The final orientation is significantly decreased as /3 is in- 
creased, e.g., for Iq ~ lO^^W/cm^, (cos0) = 0.651 and 

0. 561, and the contribution of the adiabatic ground state 
is |Cooe(0)P = 0.837 and 0.791 for /3 = 30° and 45°, re- 
spectively. For a certain angle /3, the orientation achieved 
at t = decreases as Iq is increased, cf. |Fig. 8[ Since the 
molecular dynamics of the ground state is restricted to 
the two lower pendular adiabatic states, its time depen- 
dent results for its energy and expectation values (cos^ 9) 
and (J'^) resemble the adiabatic ones. 

To illustrate the rotational dynamics of excited states, 
we show in |Fig. 9| (a) the orientation cosine (cos 9) as 
a function of l(<) for |l,l,e). The field configurations 
arc the same as in |Fig. 8| The adiabatic model pre- 
dicts a sharp wrong-way orientation. In contrast, this 
state presents a weak or even zero orientation, and in 
addition, a larger peak intensity does not imply a larger 
orientation. When the peak intensity is reached, this 
level shows a weak right- way orientation for /S — 30°: 
(cos 61) = 0.326 and 0.259, for Iq = 5 x 10" W/cm^ and 
10^^ W/cm^, respectively. For /? = 45° and the peak in- 
tensities Ig = 5 X 10^^ W/cm-^ and lO^'^W/cm^, due to 
the non-adiabatic dynamics |1, 1, e) is not oriented. 

Let us analyze in detail these results. For highly ex- 
cited states, the dynamics is more complicated. Apart 
from the doublet formation, there is another physical 
phenomenon at weak laser intensities which causes loss 
of adiabaticity. In the presence of only a weak static 
field, the M-degencracy of the states with the same 
field-free J is broken due to the quadratic Stark effect, 

1. e., the splitting goes as AE ^ Eg. As the pulse is 
switched on, the energy gap between two states of this 
./-manifold is much smaller than the rate of their cou- 
pling due to the laser field, i.e., 77 is larger than one. 
For /3 = 30°, the adiabatic parameter 77 between the 
states |l,l,e)p and |1,0, e)p, both contributing to the 
dynamics of |l,l,e), is presented in |Fig. 9] (b), and it 
achieves large values for l{t) ^ 5 x 10*^ W/cm^. As the 
states in this J-manifold are driven apart by the laser 
field, the process is non-adiabatic and there is a pop- 
ulation transfer between them. The projections of the 
time-dependent wave function in terms of the adiabatic 
states |0,0,e)p, |l,l,e)p, |l,0,e)p, and |2,2,e)p is pre- 
sented in |Fig. 9| (d) for P = 30° and Iq = lO^^W/cm^. 
Under these diabatic conditions. |Ciie(t)P decreases as 
I(t) is increased, whereas |Cioe(i)P increases. By fur- 
ther increasing I(i), the coupling between these states 
becomes very small or even zero and their energy sep- 
aration increases, so that 77 decreases and the popula- 
tion transfer is stopped. This process is so diabatic that 
the wave function does not change, but its projections 
on the adiabatic basis are modified because the basis 



9 



varies with time. For instance, the field-free state is 
|l,l,e)= cos/3|l, 1, e)p-|- sin/3|l, 0, e)p, which belongs to 
the proper irreducible representation. After swichting 
on the static field, its wave function could be approxi- 
mated by the same expression because this field is very 
weak. Once the splitting of this J-manifold is finished, 
i.e., for l{t) ~ 5 X lO^W/cm^, the contributions of the 
states |1, 1, e)p and |1, 0, e)^ are approximately sin^ /? and 

cos^ /3, respectively. Note that the states |J, A/, e)p and 
I J, M, e)p arc not related adiabatically. 
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FIG. 9. (Color online) (a) For the |l,l,e) state, expecta- 
tion value (cos 9) as a function of l{t) of a 10 ns Gaussian 
pulse. The field configurations are /3 = 30° (thick) and 45° 
(thin), with lo = 2 X 10" W/cm^ (dashed), 5 x lO" W/cm^ 
(dot-dashed) and IQ^'^W/cm^ (dotted). The adiabatic re- 
sults (solid line) are also included. For lo = 10^^ W/cm^ and 
P = 30° , (b) square of the projections of the time-dependent 
wave function on the adiabatic pendular states |l,l,e) 
(dot-dashed), |l,0,e)p (dotted), |0,0,e)p (long-dashed), and 
|2, 2,e) (short-dashed), and (c) adiabatic parameter be- 
tween the pendular states |1, 1, e)^ and |1, 0, e)^ (dot-dashed), 
|l,l,e)p and |0,0,e)p (dashed), and |1,0, e)p and |2,2,e)p 
(dotted). The dc field is fixed to Es = 300 V/cm. 

In contrast to the ground state, the wave function of 
any excited level has contributions from adiabatic states 
which correspond to different pendular doublets. As the 
laser intensity is increased, the molecular dynamics is af- 
fected by the formation of these pendular doublets. Thus, 
the final orientation could be significantly reduced com- 
pared to the parallel fields result. For instance, the time- 
dependent |l,l,e) state has contributions from the adi- 
abatic levels |1, l,e)p and |l,0,e)p, which correspond to 
the first and second pendular doublets, respectively. In 
|Fig. 9| (b) we show how the adiabatic parameters t] be- 
tween the pairs |0, 0, e)p-|l, 1, e) and |1, 0, e)p-|2, 2, e)p, 
which form the first and second doublets, respectively. 



achieve values close to 1. The final population of the 
state |l,l,e) is |Cooe(0)|2 = 0.041, |Ciie(0)|2 = 0.210, 
|Cioe(0)P = 0.603, and |C22e(0)P = 0.146, which gives 
rise to a small orientation As a consequence of this popu- 
lation redistribution to other pendular doublets, features 
of the system such as energy, alignment and hybridiza- 
tion of the angular motion do not resemble the adiabatic 
results. In particular, since the levels on the second pen- 
dular doublet possess a smaller alignment, the adiabatic 
result is larger than the time-dependent one. For in- 
stance, for /3 = 30° and Iq = 10^^ W/cm^, once the time 
evolution is finished the alignment of this state |1, 1, e) is 
(cos^ 9) = 0.879, compared to (cos^ 9) = 0.931 obtained 
for the adiabatic level |1, 1, e) . 



For /3 = 45°, despite the fact that the |1, l,e) level is 
significantly aligned, (cos^ 9) = 0.896, it is not oriented 
with (cos 61) = -0.041 for Iq = lO^^^/cm^. This state 
does not gain any orientation if stronger peak intensities 
are used. This is a consequence of the population re- 
distribution explained above. Indeed, this level could be 
considered as a dark state for the mixed-field orientation 
dynamics. This physical phenomenon is not restricted to 
this state and field configuration. We show below that 
other levels also behave as dark states. It is worth not- 
ing that if in a mixed-field orientation experiment these 
dark states form part of the molecular beam, the degree 
of orientation will be diminished. 

The population redistribution to other pendular dou- 
blets significantly affects the expectation value {cos 9s)- 
To (cos 6*3) contribute terms which mix up adiabatic 
states with different magnetic quantum numbers Since 
their wave functions could spatially overlap, their cou- 
pling matrix elements do not vanish, and (cos 9s) oscil- 
lates as t is increased. 

The phenomenon of population redistribution at weak 
laser intensities also occurs for highly excited rotational 
levels, and for them, more adiabatic states would be in- 
volved in it. Before the Gaussian pulse is turned on, the 
Stark separation of the states in a certain J-manifold is 
increased if the electrostatic field strength is enhanced. 
Then, the adiabatic parameter 77 is reduced, and the 
process of splitting of this J-manifold becomes less di- 
abatic. Indeed, for sufficiently strong dc-ficlds, the dy- 
namics would be adiabatic without population transfer 
between the states with the same field-free J. For in- 
stance, the mixed-field dynamics of the |1, 1, e) level can 
be considered as adiabatic on the weak laser field regime 
for Es > 14kV/cm and = 30°. 

Let us remark that the excited states could also suffer 
avoided crossings with adjacent levels, having different 
field-free magnetic quantum number M, as the pulse in- 
tensity is varied. The rotational dynamics in most of 
these crossings will be non-adiabatic psj . 
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A. Influence of the peak intensity Iq 

Analogously to the parallel-field configuration, we in- 
vestigate now the impact of the laser peak intensity on 
the orientation. To do so, we restrict this study to the 
following eight states: |0, 0, e), |l,0,e), |l,l,e), |2,0,e), 
|2,l,e), |2,2,e), |3,0,e), and |3,2,e). Note that they are 
related to the ones analyzed in the parallel fields configu- 
ration, by a rotation of f3 around the LFF Y axis. Their 
orientation at t = 0, i.e., (cosd) for 1(0) = Iq, is plotted 
as a function of Iq in |Fig. 10| for f3 = 30°, panels (a)-(b), 
/3 = 45°, panels (c)-(d), and /3 = 75°, panels (e)-(f), and 
Es = 300V/cm and 600 V/cni, respectively. The FWHM 
of these Gaussian pulses is fixed to r = 10 ns. 
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FIG. 10. (Color online) Expectation value {cos 6) a.t t — 
as a function of the peak intensity Iq for |0, 0, e) (red thick 
solid), |1,0, e) (gold thin solid) |l,l,e) (orange thick short- 
dashed), |2, l,e) (dark blue long-dashed) |2,0, e) (blue dot- 
short-dashed), |3,0, e) (purple dotted), |2,2,e) (green thin 
short-dashed) and |3,2,e) (pink dot-long-dashed). The field 
configurations are (a)-(b) /3 = 30°, (c)-(d) /3 = 45° and (e)- 
(f) /3 = 75°, with Es = 300 V/cm and 600 V/cm, respectively. 
The FWHMs of the Gaussian pulses are fixed to 10 ns. 

Let us start analyzing the results for the ground state. 
For all field configurations, (cos 9) shows a qualitatively 
similar behavior as a function of the peak intensity: ini- 
tially increases, reaches a maximum and decreases there- 
after. At the peak intensity where the maximum of 
{cos 9) takes place, the dynamics of this state is non- 
adiabatic, but due to the coupling of both states the ori- 
entation increases with respect to the adiabatic result. 
For a fixed peak intensity and electric field strength. 



{cos 9) decreases as /3 is increased towards 90°. For 
(3 = 75° and Iq = 2 x 10" W/cm^, the ground state 
achieves a moderate maximal orientation, (cos 9) = 0.514 
and 0.796 for Eg = 300 V/cm and 600 V/cm, respectively. 

The population transfer taking place at weak and 
strong laser intensities leaves its finger-prints in the dy- 
namics of the excited states. Compared to the paral- 
lel field results, cf |Fig. 4[ their orientation is reduced 
for any inclination angle /3 and the pendular pairs are 
not any longer formed by a right- and wrong-way ori- 
ented states. Whereas for most of the field configu- 
rations, the ground state possesses the largest orienta- 
tion, the levels |l,0,e) or |2,1,0) could achieve a sim- 
ilar or even larger orientation, e.g., for /3 ~ 45° and 
75°, lo w 5 X 10"W/cm2 and E^ = 300 V/cm. For 
f3 = 30° Eg = 600 V/cm, the degree of orientation is 
moderate for most of the states. Several dark states are 
fou nd for 1 3 = 45°: |l,l,e), |l,0,e), |3,0,e), and |3,2,e), 
see |Fig. 4| (c)-(d). For instance, the levels |l,l,e) and 
|l,0,e) are strongly aligned with {cos^ 9) = 0.927 for 
lo > 2 X lO^^W/cm^ and E^ = 600 V/cm, whereas they 
arc not orientated with {cos 9) w —0.059 and —0.043, 
respectively. For (3 = 75°, when the peak intensity of 
the Gaussian pulse is reached most of the states present 
a weak orientation, only the levels |0,0, e) and |l,l,e) 
have a large orientation for small values of lo. 

These results indicate that with a 10 ns alignment 
pulse, strong dc fields and small inclination angles are re- 
quired to reach a moderate orientation for excited states. 



B. Influence of the FWHM of the laser pulse 

For the same set of states as in the previous section, 
we analyze here how their directional properties depend 
on the FWHM of the Gaussian pulse. In |Fig. 11| (a) 
and (b) we show {cos 6) at i = as a function of r for 
respectively. The static electric field 
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lo = 10i2w/cm2 
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FIG. 11. (Color online) Expectation value {cos 9) at t = as 
a function of r. The field configurations are (a) /3 = 30° and 
(b) /3 = 45°, with lo = 10^^ W/cm^ and Es = 300 V/cm. The 
states and their labels are the same as in |Fig. lO] 

The orientation of the ground state increases approach- 
ing the adiabatic limit as t is increased, and it reaches it 
with a 50 ns pulse. We encounter several excited states 
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presenting a moderate or weak orientation, and their 
I (cos 6') I monotonically increases as r is enhanced, e.g., 
for /3 = 30° the levels |l,0,e), |l,l,e), |2,0,e), |2,l,e), 
|2,2,e) and |3,0,e) and for /3 = 45° |2,0,e), |2,l,e) and 
|2,2,e). For all of them, a 20 ns pulse is not enough 
to achieve the adiabatic regime. In contrast, other ex- 
cited levels present a very small, almost zero, orientation 
independently of the pulse duration. Some of these lev- 
els behave as dark states being strongly aligned but not 
oriented independently of the pulse duration, e.g., the 
|3,2,e) state has (cos^ 6*) = 0.755 and |(cos6')| < 0.04 for 
/3 = 30° and any value of r. An analogous behavior is 
found for the levels |1,0, e), |l,l,e), |3,0, e) and |3,2,e) 
at /? = 45°. As described above, this phenomenon is due 
to the non-adiabatic dynamics at weak laser intensities 
when the levels of the J- manifold are driven apart, and 
it takes places even for 50 ns pulses. 



C. Influence of the electrostatic field strength 

In this section, we consider two inclination angles, and 
investigate the impact of the electrostatic field on the 
mixed-field orientation dynamics of the same states. Fig- 
ures [T2] (a) and (b) illustrate the behavior of {cos 9) at 
f = as a function of Eg for /3 — 30° and 45°, respectively. 
The laser pulse has t = 10 ns and lo = 10^^ W/cm^. 
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FIG. 12. (Color online) Expectation value of {cos 9) at t = 
as a function of Eg for the field configurations r = 10 ns, 
lo = 10^^W/cm^ and (a) /3 = 3 0° and (b) /? = 45°. The 
labeling of the states is done as in |Fig. 10| 



The ground state presents the largest orientation, 
which increases as Eg is enhanced, being strongly ori- 
ented for sufficiently large fields, e.g., {cos9) > 0.9 for 
Es > 600V/cm and /3 = 30°. Regarding the excited 
states, their orientation strongly depends on the incli- 
nation angle. For /3 = 30°, |(cos6')| monotonically in- 
creases till it reaches a plateau-like behavior, and they 
show a moderate orientation. Indeed, for (3 = 30°, 
Ig = lO^^W/cm^ and Eg — 2kV/cm, we obtain at the 
maximum of the Gaussian pulse {cos 9) — —0.402 and 
0.365 for the states |1, 0, e) and |1, 1, e), respectively. For 
P = 45°, the level |2,l,e) presents a large orientation: 
(cos 61) > 0.8 for Eg > 800V/cm. There are some darks 
states for /3 = 45°, which are not oriented even when dc 
fields of 2 kV/cm are used, e. g., |3, 0, e) and |3, 2, e). 



For non-parallel fields, a strong dc field does not ensure 
a large orientation for excited rotational states. If the 
aim is a strongly oriented molecular ensemble, then this 
should be as pure as possible in the ground state. 

In the Hamiltonian |(3)[ the term — /xEg sin /3 sin cose/) 
is responsible for the mixing of states with different field- 
free magnetic quantum numbers. On the weak dc field 
regime, the mixing between these states is so small that 
M could be considered as conserved, and this term could 
be neglected. By increasing Eg, this coupling between 
levels with different field- free M becomes important, and 
this should affect the molecular dynamics. Thus, the 
questions that remain open is how important is the X- 
component of the electrostatic field to the dynamics, and 
for which electric field regime, we could only consider its 
Z-componcnt Eg = Eg cos /3Z. 

As indicated above, even for tilted fields, the dynamics 
of the ground state can be described by a two state model. 
Its energy separation to the next state with M 7^ is 
of the order of B and larger than the dc field coupling 
to these levels. Thus, for Eg < 20kV/cm, the dynamics 
considering the dc field is equal to the one obtained when 
only its Z-component is included. 

For the excited states, the answer to these questions 
depends on how the initial wave function, before the 
fields are switched on, is constructed. The first op- 
tion is to proceed as indicated at the beginning of this 
section; the field-free /3 7^ 0° and /3 = 0° wave func- 
tions are related by a rotation of /3 around the Y-axis 
\J,M,e)=RY{l3)\J,M,e)°. In this case, for the level 
|l,l,e), some differences in its orientation are observed 
for Eg > IkV/cm with Kcos^)! being larger if the two 
components of Eg are considered. These differences are 
augmented as Eg is increased, e. g., for a 10 ns laser pulse 
with lo = 1 X 10^2 W/cm^, p = 45° and Eg = 5kV/cm, 
we obtain at f = (cos 6') = —0.629 compared to 
{cos9) = —0.019 if only the Z-component of Eg is in- 
cluded. By increasing Eg this state will achieve an adia- 
batic dynamics only if both components of the static field 
are present. The second option is to construct the field- 
free /? 7^ 0° wave function equal to the field-free /3 = 0° 
one. In this case, the results resemble those of the paral- 
lel field configuration taking into account cos p as scaling 
factor for the static field strength. 



D. Influence of the inclination of the fields 



The symmetries of the rotational Hamiltonian (1) 



see 

ISec. Ill and, therefore, the rotational dynamics strongly 
depend on the angle between the fields. In this section, 
we investigate in detail the impact of the inclination angle 
in the mixed-field orientation dynamics. 

For the ground state |0, 0,e), the orientation cosines 
{cos 9) and (cosflg) are plotted in |Fig. 13[ func- 
tion of /3, together with the adiabatic results. For a 
weak dc field and strong laser field, the following rela- 
tion (cos6'g) « (cos 61) cos /3 is satisfied within the adi- 
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abatic limit. In (cos 0s) the term (sin cos c/i) sin /3 has 
been neglected, which can be done as far as the mix- 
ing between states with different field-free M is very 
small. By increasing the electrostatic field strength, a 
regime would be encountered where this approximation 
does not hold any longer. An analogous relation is sat- 
isfied between the time-dependent orientation cosines of 
the ground state. For Iq = 2 x 10^^ W/cm^, its orienta- 
tion {cos 9) shows a plateau-like behavior till (3 = 50°, 
which is very close to the adiabatic limit. By further in- 
creasing /?, {cos 9) decreases and approaches to zero. For 
/3 = 90°, the states in a pendular doublet have differ- 
ent symmetry and are not coupled by the dc field, thus 
they might be strongly aligned but not oriented. For 
lo = lO^^W/cm^, {cos 9) monotonically decreases as (3 
is increased towards 90°, and its value is always smaller 
than for Iq = 2 x lO^^W/cm^. For both laser fields, 
(cos 9s) decreases as /3 is increased. 




30 45 60 

PC) 



FIG. 13. (Color online) Expectation values {cos 6) (thick solid 
line) and (cos^s) (thin sohd hne) at t = as a function of 
P for the ground state. The peak intensities are (a) Iq — 
2 X 10" W/cm^ and (b) Iq = 10^^ W/cm^ The adiabatic 
results for {cos 9) (thick dashed line) and (cosSs) (thin dashed 
line) are also presented. The FWHM of the laser pulse is fixed 
to r = 10 ns and the dc field to Ea = 300 V/cm. 

In Figs. [13] (a), (b), (c) and (d), we present the orien- 
tation cosine of the pairs |0, 0, e)-|l, 0, e), |1, 1, e)-|2, 1, e), 
|2, 2, e)-|3, 2, e) and |2, 0, e)-|3, 0, e), respectively, as a 
function /3. The static field strength is Eg = 300 V/cm 
and we consider two Gaussian pulses of 10 ns FWHM and 
peak intensities Iq = 2 x 10^-^ W/cm^ and lO^^W/cm^. 
Due to the complicated field-dressed dynamics of excited 
states for 0° < /3 < 90° with contributions from sev- 
eral pendular pairs, in (cos 6*5) the term (sin 6* cos 0) can- 
not be neglected. Then, the simple relation (cos 0s) ~ 
(cos 9) cos f3 does not hold for these levels. 

Based on the adiabatic theory, the ground state and 
the level |1, 0, e), should present the same orientation but 
with opposite directions. However, this is only satisfied 
for (3 = 0°. Due to the non-adiabatic effects at weak 
laser intensities, its | (cos 0)1 is smaller than the corre- 
sponding value of |0, 0,e) for < f3 < 90°. For the 
second doublet, cf. |Fig. 14| (b), {cos 9) oscillates as /3 
is varied, and the orientation even changes its direction. 
Both states could present a moderate orientation at a 
certain value of /?. The pendular regime is not achieved 
by the third and fourth pairs with a 10 ns laser pulse and 
Iq = 2 X 10^^ W/cm^. and their orientation is either zero 




FIG. 14. (Color online) Expectation value (cos 6) a.t t = 
as a function of /3 for the states (a) |0,0, e) (black), and 
|l,0,e) (blue), (b) |1, l,e) (black), and |2, l,e) (blue), (c) 
|2,0,e) (black), and |3,0,e) (blue), (d) |2,2,e) (black), and 
|3, 2,e) (blue). The field configuration is Es = 300 V/cm and 
lo = 2 X 10" W/cm^ (solid lines) and 10^^ W/cm^ (dashed 
lines). The FWHM of the laser pulse is fixed to r = 10 ns. 



or very small independently of f3. For Iq = 10^^ W/cm^ 
and (3 = 0°, these four states show a moderate orienta- 
tion, which is reduced for any other angle, being small for 
(3 > 60°. At the strong peak intensity Iq = 10^^ W/cm^, 
in all pendular doublets one of the two levels presents the 
dark behavior with respect to the mixed-field orientation 
dynamics at a certain angle (3. 

These results show that if the molecular beam is rota- 
tionally cold, a small inclination angle will optimize the 
degree of orientation observed in the experiment. 



VI. ROTATIONAL DYNAMICS ONCE THE 
LASER PULSE IS TURNED ON 

Let us investigate the dynamics for t > assuming 
that the laser peak intensity, reached at t = 0, and the 
dc field strength are kept constant for t > 0; i. e., l{t) = Iq 
and Es(i) = Es for t > 0. At t = 0, the time-dependent 
wave function can be expressed in terms of the corre- 
sponding adiabatic basis. Since the Hamiltonian is time- 
independent for t > 0, the contribution of each adiabatic 
state remains constant as t is increased. For a certain 
state I7), the expectation value of an operator A in this 
adiabatic basis reads as 



+ 2^|Q(0)||Cfe(0)|p(7,|i|7fe)pCos 
j<k 



(7) 



with Cj{0) being the weight at t = of the adiabatic 
state |7j)p to the wave function of I7), AEjk the energy 
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splitting between the adiabatic levels |7j)p and |7fc)p, and 
djk the phase difference of Cj{0) and Cfe(O). 

Based on the results presented above, the time- 
dependent wave function could have contributions from: 
i) only the adiabatic levels forming a pendular doublet, 
or ii) several adiabatic levels from at least two pendu- 
lar doublets. All the states for (3 = 0°, and the ground 
state for 0° < /3 < 90° could belong to the first case. 
Whereas, the second one refers to all excited states when 
0° < (3 < 90°, unless the static field is very strong. 
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FIG. 15. (Color online) Orientation cosines once the peak 
intensity and dc field strength are kept constant for t > 0. 
For the state |l,l,e), (a) expectation value {cos 9} with 
lo = 2xlO"W/cm2 and /3 = 0° (solid), 30° (dashed) 
and 45° (dotted); and (b) (cose^) for /? = 30° and 
lo = 2 X 10" W/cm^ (solid), 5 x lO^W/cm^ (dashed) and 
10^2 W/cm^ (dotted). The dc field is = 300 V/cm. 

Let us first analyze the case when the dynamics takes 
place within a pendular doublet. If the adiabatic states 
arc not fully oriented, the coupling term in Eq. (7) is 



non-zero and this expectation value oscillates for t > 
with the frequency equal to the energy splitting of the 
corresponding pendular doublet. For the |l,l,e) state, 
this behavior is shown for the time evolution of (cos 6) in 
|Fig. 15| (a), with Iq = 2 x 10"W/cm2, = 300 V/cm 
and /3 = 0°. An analogous behavior is obtained for the 
ground state and 0° < /3 < 90°. By further increasing 
the peak intensity, the orientation of the adiabatic states 
increases, the coupling terms are reduced approaching 
zero, and these regular oscillations will disappear. 

When two pendular doublets participate in the dynam- 
ics, this oscillatory behavior becomes irregular, because 
the frequencies associated with the energy separations 
within each pendular doublet and between two of them 
do not form a commensurable set. As an example, we 
show in |Fig. 15| (a) these irregular oscillations of (cos 9) 
for |l,l,e) with /3 = 30° and 45°, E^ = 300 V/cm and 
lo = 2 X 10^^ W/cm^. By increasing Iq, the dynamics 
of this state still has contributions from different pendu- 
lar doublets, but the two states in a pendular pair are 
not populated. As a consequence, the coupling terms are 
reduced and the oscillation decreases or even disappears. 

For 0° < /3 < 90°, the time evolution of (cos 9s) is dom- 
inated by the couplings of adiabatic levels from doublets 
with \ AM\ sa 1. This is illustrated in |Fig. 15| (b), for the 
state |1, 1, e), Iq = 2 x 10" W/cm^, 5 x 10" W/cm^ and 
10^2 W/cm^ and /3 = 30°. Independently of Iq, in this 
time scale (cos6's) oscillates with the largest frequency 



given by the energy gap between of the two pendular 
doublets involved, which is similar for the three peak in- 
tensities. On a larger time scale, the frequencies due to 
the states in a doublet will modulate the oscillations of 
(cos 9s) in the weak field regime. 



VII. SWITCHING ON THE LASER PULSE 
FIRST: ORIENTATION OF THE ALIGNED 
PENDULAR STATES. 

In previous sections, the field configuration was based 
on the mixed-field orientation experiments [13, [13, [1^ ■ 
Here, we investigate the molecular dynamics when the 
temporal order of the fields is inverted: the Gaussian 
pulse is switched on first, its peak intensity is kept con- 
stant, and then the static electric field is turned on. 
While the laser field is switched on, the pendular dou- 
blets of quasi-degenerate states with opposite parity are 
formed. This process is adiabatic [4l|, [4l|, and for a suffi- 
ciently large peak intensity, these two levels are strongly 
aligned but not oriented. By turning on the static field, 
these states have the same symmetry and they should be 
oriented due to their interaction with this field. For this 
field configuration, we check now the validity of the adi- 
abatic predictions (Tsl . [l9| by comparing them to a time 
dependent analysis. 

The peak intensity of the Gaussian pulse is reached 
at < = and kept constant afterwards. At this point, 
if Iq is large enough, the energy gap between the states 
in a pendular doublet is much smaller than the energy 
gap with the neighboring doublet. Then, for a certain 
pendular level, its rotational dynamics for i > 0, i.e., 
when the static field is switched on, could be approx- 
imated by a two-state model involving the two levels 
forming the corresponding pendular doublet [l9| . At 
t = 0, i.e., Es(0) = and 1(0) = Iq, the pendular states 
are \ipi) with I = e and o indicating even or odd par- 
ity. Under this approximation, we assume that the levels 
Hr) = ^ i\ipe) + iV'o)) and Itp^j) = (He) " \ipo)) are 
right- and wrong-way oriented, respectively. The two- 
state-model Hamiltonian yields as 



H(t) 





—fiVst{cos t 



-fJ,Vst{cOS 9s 

AE 



where we have taken Es{t) = Vgt with Vs = Eg/To and Tq 
being the switching on speed and time, respectively. This 
time To is chosen so that if these states are exposed only 
to this field, the turning-on process is adiabatic. Note 
that wc have taken (i/>e|^f|V'e) = 0, {^o\H\^po) ~ AE, 
and {ipe\H\ipo) = {iJo\H\ipe) = -fiVst{ipe\cos9s\ipo) = 
-~liVst{cos6s)eo- The time-dependent Schrodingcr equa- 
tion associated to this Hamiltonian admits a scaling fac- 
tor. That is, when the dynamics is adiabatic using Vs for 
a pendular doublet with energy splitting AE at 1(0) = lo, 
then, for a peak intensity Iq and AE' = kAE, the dy- 
namics is adiabatic for v' = k'^Vs- 
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FIG. 16. (Color online) For the ground state, (a) expectation 
value {cos 8} and (b) weight of the adiabatic ground state 
on its time-dependent wave function as a function of Es(t), 
for turning on speeds Vs = 10^° (orange dashed), lO'' (red 
solid), 10* (blue dotted) and lO'^ V/scm (pink dot-dashed), 
and adiabatic results (thin solid) The fields are parallel and 
the Gaussian pulse has r = 10 ns and lo = 2 x 10^^ W/cm^. 



For the sake of simplicity, we focus on the ground state 
in a parallel-field configuration. For several switching on 
speeds, |Fig. 16| (a) and (b) display the directional cosine 
and the population of the adiabatic ground state, respec- 
tively, as a function of Es(t). The Gaussian pulse has 
10 ns FWHM and peak intensity lo = 2 x 10" W/cm^. 
Before the dc field is turned on, the alignment of the 
ground state is (cos^ 9) = 0.845, the energy separation 
within this pendular pair is AE « 5.36 x lO^"' cm^^, 
and there are 2.4 cm~^ to the next pendular doublet. 
For Eg ~ IV/cm, the coupling term is fi{cos9)eo^s = 



1.09 x 10"'=' cm"^ with (cos( 



0.915. In an adiabatic 



picture, the energy gap AE' can not be neglected, and 
as Es{t) is increased the ener gy o f the ground state does 
not increase linearly with Eg [l^. For Vs = 10^° V/scm, 
the adiabatic parameter rj 1, the rotational dynam- 
ics is non-adiabatic and there is a population transfer 
between the two states in this doublet. We note that 
for this process, the adiabatic parameter ry is defined 
as in Eq. (6) but replacing the laser field interaction 
^L(t) 1(4)1 by the dipole term Hs{t) [(3)] The ground 



state presents a moderate orientation bellow the adi- 
abatic limit due to the contributions of the adiabatic 
states |0,0, e)p and |1,0, e)p, |Cooe(i)P decreases until a 
minimum value showing a smooth oscillation afterwards, 
cf. |Fig. 16| (b). Due to the coupling term, {cos 9) os- 
cillates as Es(t) is increased, and its frequency is equal 
to the energy separation between the adiabatic levels 
|0, 0, e)p and |1, 0, e)p. A similar behavior is observed for 
10^ V/scm, but the orientation of the ground state oscil- 
lates around a value closer to the adiabatic limit because 
the process is more adiabatic and |Cooe(i)P ~ 0.956 for 
Es(t) > 100 V/cm. For = 10^ and lO'^ V/scm, the 
dynamics can be considered as adiabatic, with |Cooe(i)P 
being larger than 0.999. However, for Vg = 10^, {cos9) 
still oscillates around the adiabatic value. 

By increasing the peak intensity of the laser pulse, the 
energy splitting of the levels in a pendular doublet is de- 
creased, but their coupling due to the dc field is not signif- 
icantly modified. Thus, the rotational dynamics becomes 
more diabatic, and larger turning-on times are needed to 



achieve the adiabatic limit. For Iq = 5 x 10"'^^ W/cm^, 
the ground state is separated by AE w 7.7 x 10~^ cm~^ 
from |1,0, e)p, and by 3.9 cm~^ from the next pen- 
dular doublet. The coupling due to the dc field is 
^{cos9)eoEs = 1.13 X 10"^ cm-1 for E^ = IV/cm and 
with (cos 6')eo = 0.948. According to the scaling law of 
the time-dependent Schrodinger equation, this process 
would be adiabatic for a speed of Vg < 20 V/scm. 

For lo = 10^^ W/cm^ and the ground state, wc find 
AE « 3.3 X 10"^° cm-\ 5.7 cm"! to the second dou- 



blet, and /x(cos( 



,E, = 1.18 X 10" 



-1 for E, = 



IV/cm and with (cos 6')eo = 0.964. Within an adia- 
batic framework, as Eg is increased the ground state 
energy can be approximated by the pseudo-first-order 
Stark linear effect I^E^ p (0, 0, e|cos6'|l, 0, e) | Note 
that AE is smaller than the dc filed coupling even for 
Eg ~ 10^"* V/cm. Based on the scaling law of the two- 
state model Schrodinger equation, the dc field should be 
turned on very slowly, Vg ^ 10~^ V/scm, to achieve the 
adiabatic limit. For larger turning-on speeds, the dy- 
namics is so diabatic that the |0, 0, e) wave function does 
not change, and its projections on the adiabatic states 
|0,0, e)p and |1,0, e)p are close to the field-free values 
even for Vg ~ 10'"' V/scm. 



VIII. CONCLUSIONS 

In this work, we have investigated the mixed-field ori- 
entation dynamics of linear molecules. The richness and 
variety of the field-dressed rotational dynamics has been 
illustrated by analyzing in detail the directional proper- 
ties of several low- lying states. In particular, we have 
explored the degree of orientation as the peak intensity 
and FWHM of the Gaussian pulse, the electrostatic field 
strength and the angle between both fields are varied. 

By considering prototypical field configurations used 
in current mixed-field orientation experiments, we have 
proven that the assumption of a fully adiabatic dynam- 
ics is incorrect. For parallel fields, a non-adiabatic trans- 
fer of population takes place when the quasi-degenerated 
pendular doublets are formed as the laser intensity is 
increased. As a consequence, the time-dependent re- 
sults for the degree of orientation arc smaller than the 
predictions of the adiabatic theory. Using current avail- 
able experimental peak intensities, longer laser pulses or 
stronger static fields will increase the degree of orienta- 
tion even for highly excited states. In particular, we have 
provided the field parameters under which the mixed- 
field orientation dynamics will be fully adiabatic. We 
have also shown that the field-dressed dynamics is more 
complicated if both fields arc tilted. Apart from the non- 
adiabatic effects when the pendular doublets are formed, 
at weak laser intensities there is also population trans- 
fer due to the splitting of the states within a J-manifold 
having now the same symmetry. For non-parallel fields, 
we have shown that the ground state is strongly oriented, 
whereas excited states might only present a moderate ori- 
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entation, and, furthermore, some of them could behave as 
dark states to the mixed-field orientation dynamics. The 
requirements for an adiabatic dynamics arc now more 
difficult to satisfy for excited levels than for the ground 
state. Again, we have indicated the field configuration 
that will give rise to an adiabatic mixed-field-orientation. 
If the peak intensity is kept constant after turning on the 
pulse, we have shown that the orientation of the states 
might oscillate with time due to the non-adiabatic dy- 
namics. Finally, we have investigated the molecular dy- 
namics when the temporal order of the fields is inverted. 
We have shown that once the ground state is adiabati- 
cally aligned, the switching on of the dc field has to be 
very slow to achieve a significant orientation. 

Although our study is restricted to the OCS molecule, 
we stress that the above-observed physical phenomena 
are expected to occur in many other polar molecules. In- 



deed, the Hamiltonian can be rescaled, and the above 
results used to describe another polar linear molecule. 
In addition, due the complexity of their rotational level 
structure of asymmetric tops, these non-adiabatic effects 
should have a negative impact in their mixed-field orien- 
tation experiments p^ . [25| . 
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